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Below pion threshold virtual Compton scattering off the nucleon gives access to 
the generalized electromagnetic polarizabilities. Different theoretical results for 
the generalized polarizabilities have been compared, fn particular, the influence 
of the generalized polarizabilities on the unpolarized cross section and on double 
polarization observables has been investigated. Predictions for these observables 
have been obtained in the linear sigma model and chiral perturbation theory. 

1 Introduction 

Compton scattering has a long history as an interesting tool to investigate the 
nucleon’s structure. Though it has been applied in different kinematical re¬ 
gions, Compton scattering is of particular interest below the threshold of pion 
production, where the electromagnetic polarizabilities of the nucleon can be 
studied. While in the past the main focus was on real Compton scattering 
(RCS), recently much effort has been devoted to virtual Compton scattering 
(VCS), 7* -b A —> 7 -b A, with a virtual photon 7* in the initial state and a 
real photon 7 in the final state. In comparison to RCS, the VCS reaction con¬ 
tains much more information because of the variable mass and the additional 
longitudinal polarization of the virtual photon. 

As shown in Fig. ^ VCS can be realized experimentally by electron-nucleon 
bremsstrahlung, e-b A ^ e'-b A'-by. Besides the VCS mechanism, the brems- 
strahlung contains the Bethe-Heitler scattering (BH), where the real photon 
is emitted from the incoming or outgoing electron. Therefore, the total am¬ 
plitude T®® ^ is the coherent sum of the VCS amplitude gjj 
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Figure 1: Feynman diagrams of the reaction e{k) + N{p) e'(k') + N'(p') + 7(g0- (^)- 
VCS contribution containing the Born and residual amplitudes, (b) and (c): Bethe-Heitler 
scattering. 


amplitude 

rpee'^ rpVCS _j_ rjiBH 

If not stated differently, we use the cm frame of the final state photon and 
nucleon throughout this article. The 4 -momenta of the incoming photon and 
nucleon are denoted by = {lu, q) and = {E, —q), respectively, while the 
outgoing particles are characterized by q'^ = (w', q') and p'^ = {E', —q'). We 
also define = —q^ and q= |(f|. For the following discussion it is convenient 
to introduce the quantity 


Qo = = 2mAr 



( 2 ) 


According to Eq. (|^, the variable Qq can be replaced by q and vice versa. 

Below pion threshold, the analysis of both RCS and VCS is based on low 
energy theorems, which rely on fundamental symmetries as Lorentz invariance, 
gauge invariance and crossing symmetry. The decomposition of the scattering 
amplitude into a Born contribution and the residual part plays an important 
role in the derivation of these theorems. In the case of VCS we have 


rjnvcs _ j'VCS _|_ xvcs ^ 3 ^ 

where the Born amplitude shows the typical l/ta' singularity of brems- 

strahlung, wh^ the residual amplitude is proportional to co' in the low 

energy limit EJjj. The Born amplitude is entirely determined by the nucleon 
mass, its charge and the electromagnetic form factors, Ge and Gm, he., by 
definition contains only properties of the nucleon in its ground state. 

The influence of the excitation spectrum enters through the residual ampli¬ 
tude which mn be parametrized by the generalized polarizabilities (CPs) of 
the nucleon El. To leading order in w', there appear 10 GPs in These 

CPs depend on Qq or q and are related with the electromagnetic polarizabil¬ 
ities of RCS whose leading contribution is completely determined by the well 
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known electric (a) and magnetic {(3) polarizabilities. We discuss the GPs in 
the following section in more details. 

The measurement of the GPs requires data with a very high accuracy 
which can only be obtained by means of the new electron accelerator facilities 
like MAMI, MIT-Bgdes and Jefferson Lab. The VGS reaction has already been 
measured at MAMlH for q = 600 MeV, it will be investigated at g = 240 MeV 
at MIT-Bates □, while the activities at Jefferson LabB will concentrate on the 
region of g « 1 GeV. Very recently, the first prelimin^y results for a(g) and 
/3(g) have been extracted from the MAMI experiments. 

On the theoretical side the GPs were predicted by various approaches. Calcu¬ 
lations of all GEa have been performed in the nonrelativistic constituent quark 
model (CQM) ElEl, chiral perturbation theory in the heavy baryon foriri^an 
(HBChPT) including the A resonance as a dynamical degree of freedopifl’EHl3, 
and in the one-loop approximation of the linear sigma model (LSM)lEl’t3. The 
g behaviour of a and /3 was determined in an effective Lagrangian model 113 
containing several resonances and the exchange of tt^ and tr mesonSpin the t 
channel. In addition, there exists a Skyrme model prediction of a(g)li3, while 
the paramagnetic part of /3(g) was emulated in a relativistic quark-model for¬ 
mulated in the light-front dynamics 113 . 

The influence of the GPs on observables as the unpolarized cross section of the 
reaction p{e,e'p')^ and asymmetries for beam-recoil polarization, p(e, 
has been studied in the effective Lagrangian approach by Vanderhaeghenll3ii3. 
In these studies the contributions to the amplitude were investigated to 

all orders in oj', i.e., the effects of higher order GPs were automatically incor¬ 
porated. In contrast to this, we will only consider the dipole approximation 
to . Regarding the validity of this approximation, see the discussion in 

sec. 2.1. In the unpolarized case, the w' dependence of the effect due to the 
GPs was previously calculated in HBGhPTllll. In addition to unpolarized cross 
sections, we also predict asymmetries of beam-target polarization, p(e, e'p' )j, 
which can be explored at Bates using the BLAST target-detector system. We 
will show predictions of the LSM and HBGhPT focusing on the influence of 
the GPs as function of the 7*7 scattering angle. 


2 Generalized Polarizabilities 

2.1 Definitions and Constraints 

The concept of GPs was first introduced in connection with nuclear targetsU 
and has recently been discussed in detail for the specific case of the nucleonEl. 
The GPs can be defined through the multipoles ^ g) of the resid¬ 

ual amplitude. In the notation of the multipoles, p {p') denotes the type of the 
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Table 1: Possible transitions to VCS in the case of the dipole approximation for the real 
photon (E: electric, M: magnetic, C: Coulomb). The last column gives the notation of the 
GPs. 


final photon 

initial photon 

S 

GPs 

El 

Cl 

0,1 

p(01,01)S’ 

El 

El 

0,1 

p(01,01)S p(01.1)S 

El 

M2 

1 

p(01.12)l 

Ml 

Ml 

0,1 

p(ll,ll)S’ 

Ml 

CO 

1 

p(11.00)l 

Ml 

C2 

1 

p(11.02)l 

Ml 

E2 

1 

p(11.02)l p(11.2)l 


incoming (outgoing) photon (p = 0: Coulomb, p = 1: magnetic, p = 2: elec¬ 
tric) while L {L') refers to the angular momentum of the initial (final) photon. 
The quantum number S indicates a no spin-flip {S = 0) or a spin-flip {S = 1) 
transition. 

In order to define the GPs it is suitable to replace all electric transitions. As a 
consequence one introduces the socalled mixed multipoles ^ describ¬ 
ing a well defined mixture of an electric and a charge transition in tljije-|initial 
state. More details on this technique can be found in the literature S’EJ. The 
GPs are finally given by 


pip'L',pL)S_ 




^(p'L'.pL)5(^,^_) 


Lj'—O 


^iL'qL+i 




(p, p' = 0,1) 

ip' = 0,1), 


( 4 ) 


uj'—O 


and carry the same quantum numbers as the corresponding multipolesB. 

In the dipole approximation for the real photon [L' = 1), selection rules due to 
parity and ang ular momentum conservation lead to the ten possible transitions 
listed in Table 3 scalar GPs {S = 0) and 7 vector GPs {S = 1). Note that 
the dipole approximation is equivalent to keeping only the leading (linear) 
order in uj' for the residual amplitude. Since the relevant expansion parameter 
is io'/rriTr-, the dipole approximation describes the experimental data for RGS 
only if ui' < 60 MeV, well below the pion mass m^n-. However, in VCS the 
situation is more complicated. As has been pointed out recently, higher order 
teiHis in iu' become more important if q is of the same order of magnitude as 
w'Ej. Therefore, beyond the restriction given by the pion mass, the condition 
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oj' should be fulfilled to guarantee the validity of the dipole approximation 
in VCS. 

In part the GPs are related to the electromagnetic polarizabilities which govern 
the low-energy expansion of the RCS amplitude. The scalar polarizabilities 
p(oi’Oi) 0 (g) and generalize a and /3 to the case of virtual photons, 


with e^/Air « 1/137. Two of the vector GPs ar^onnected with the four spin 
polarizabilities 7 ^ of RGS as defined by Ragusac2l. The relations readcil 


The remaining two combinations of spin polarizabilities of RCS can not be 
related to the GPs in the dipole approximation. In particular, the forward 
spin-polarizability 7 = 71 — 72 — 274 is not contained in this kinematical limit 
ofVCSEil. 

As has been shown by ouK-explicit calculation in the LSMIIj and later by 
a model-independent analysis 123 based on charge conjugation symmetry and 
nucleon crossing, the ten GPs have to satisfy four constraints. For instance, 
the three scalar GPs obey the condition 


/_p(01.1)0 

47r 


( 7 ) 


3mNq^ 


a{q) + P{q) 


(7) 


Similarly, there are three relations among the vector GPsEj’Ej. In addition 
to these relations, the vector GPs fulfill specific constraints at g = 0. Three 
vector GPs vanish at the origin, 


p(oi.oi)i(o) = p(ii4i)i(o) = p(ii.oo)i(o) ^ 0 ^ (8) 


while the remaining four satisfy the condition 

p(oi. 12 ) 1 ( 0 ) -f - 2^5171 ArP(“’^)^(0) = 0. (9) 

Because of the constraints, only six instead of ten GPs are independent. Though 
there is some arbitrariness in the selection of the independent GPs, a natural 
choice h^ been proposed in a recent review article of Guichon and Vander- 
haeghenES. The authors eliminate the three mixed GPs p(oi’i)o ^ p(oi.i)i ^ 
p(ii, 2 )i quantity In the following discussion we focus on the 

remaining six GPs. 
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2.2 Results and Discussion 


Both nonresonant and resonant excitations of the nucleon contribute to the 
GPs. However, in the kinematical region of Compton scattering below pion 
threshold, the nonresonant s-wave production of pions is of particular impor¬ 
tance. This process is determined by chiral symmetry and can be described by 
pion-loop diagrams in ChPT or in the LSM. On the other side, we also expect 
large influences of resonances. In particular, the A(1232) and the Di3(1520) 
resonances, which are both clearly visible in the photoabsprption spectrum of 
the nucleon, signihcantly contribute to some of the GPsBEA. 

In addition to these excitations of the nucleon, the chiral anomaly (tto exchange 
in the t channel) is quite essential in the case of the vector GPs. At low val¬ 
ues of q, the shape for three of the four independent vector GPs is strongly 
governed by the anomaly El’tHi, as will be discussed in sect. 3. 

Before comparing the numerical predictions of different apprnarbps we 
discuss the analytical results for th^GPs obtained in the LSMO’EJ and in 
HBChPT to lowest order (0(p^)) HEEI. The chiral expansions of the GPs at 
Qq = 0 contain the most instructive information. As an example, we find 
for the electric and magnetic polarizabilities of the proton in the case of the 
one-loop approximation to the LSM 


ap( 0 ) 




Stt 


1927r3m^ i2fj, 
13.6- 8.8-f 4.2 


18 In/j, - 

+ Oifj.) 


33 

T 


Oifi) 


X IQ-'^ fm' 


/3p(0) 




^ -I- 181n/r -I- ^ -I- 0{fi) 


192TT^m% lifj, 

1.4-8.8-k8.1-f C)(m) 


X lO""* fm' 


( 10 ) 


with fi = mT^/iriN and = TaNgA /~ 13.4 the pseudoscalar pion-nucleon 


coupling constant. It is interesting to note that the chin 
agree with a one-loop calculation of relativistic GhPTC 


expansions of Eq. (|10D 
to all orders in g. The 


corresponding results of HBChPT to 0(p^) are given by the leading term in 

phtained for all of the GPs 


Eq. (|^). Similar chiral expansions have b^ 
and for their derivatives with respect to Qq Both chiral approaches 

completely agree in the sense that in all cases the leading contribution of the 
expansion in the LSM is exactly the result of HBChPT to 0{p^). 

The chiral expansions of the polarizabilities show a common feature: compared 
to the leading term the contribution next to leading order has a different sign. 
Such an alternating behaviour of the chiral series can also be observed in 
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HBChPT to 0{p‘^) for q;(0) and /3(0) i^ne neglects counterterm contributions 
which arise from resonance saturation O. 

The pion cloud of the nucleon causes a diamagnetic response and therefore a 
negative magnetic polarizability. According to Eq. this diamagnetic part 
of (3 is not contained in the leading chiral contribution but in the terms next 
to leading order (logarithmic plus constat term). In Fig. ^we plot th^(^s as 
a function calculated in the CQM0, in HBChPT to order and 

in the LSMEjEj. We first consider the scalar polarizabilities a and (3. At the 
real photon point (Qq = 0), the LSM underestimates the empirical values 
(a“P(0) = 12.1 ± 1.0 X 10-4 fm^, /3 “p( 0) = 2.1 =F 1-0 x 10-4 fm^) to some 
degree. This drawback can be attributed to the neglect of t-channel exchange 
of heavier or more meson states, and of nucleon resonances in the s-channel, 
e.g. Di 3(1520) in the case of a and A(1232) with its strong paramagnetic 
(quark spin-flip) contribution to (3. 

In HBChPT the results at the real photon point are in good agreement with the 
experiment. However, at least for the magnetic polarizability, this agreement 
is somewhat accidental because both the paramagnetic contribution of the 
A(1232) and the main diamagnetic contribution of the pion cloud have not been 
taken into account. The results for a and 13 to show a large cancellation 

between diamagnetic and paramagnetic terms Ej, where the final results for 
the proton are ap = 10.5 ± 2.0 x 10-4 fm^ and f3p = 3.5 ± 3.6 x 10-4 fm^. In a 
recent development the A was incorporated in to the effective chiral lagrangian. 
Unfortunately, this formalism predicts too large numbers for a and /3 in a 
calculation to C>(e^)0E3. 

Though the values of a and (3 at the real photon point are different for HBChPT 
and the LSM, both chiral calculations predict a very similar behaviour for 
these polarizabilities. In the CQM the Qq dependence is rather different in 
comparison to the chiral calculations, the discrepancy being strongest in the 
case of f3. 

At low Qg, the polarizability is very large in the CQM in contrast 

to the chiral calculations. This behaviour is due to the neglect of diamagnetic 
contributions in the CQM. Furthermore, the CQM predicts finite values for 
the CPs origin and therefore violates the model- 

independent constraints of Eq. (^. The most remarkable discrepancy between 
the CQM and the chiral predictions is in the absolute values of 
p^ii, 02 )i (jjffgp ]gy ^.y^ro ordors of magnitude. 


3 Observables 
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Figure 2: Comparison of GPs as function of Qq obtained in different models. The anomaly 
is not included. Solid line: LSM, dashed line: CQM, dash-dotted line: HBChPT to 0{p^). 
Note that the CQM results have been scaled. 









3.1 Formalism 


In the case of the unpolarized reaction p(e, e'p' )7 we consider the five-fold 
differential cross section, 


(f’C 


^labid^p’ )c 


= T. 


ee^7 |2 


( 11 ) 


where Ki represents a phase space factor. The cross section depends on the 
five variables 

q, 0, <t), e, (12) 

with 9 the angle between the two photons, (j) the azimuthal angle between the 
leptonic plane and the reaction plane, and e the transverse polarization of the 
virtual photon. According to Eqs. (0) and (i, the scattering amplitude T®® 
consists of three parts, 

J.ee'7 ^ r^BH j.VCS j^VCS ^ 

As explained in sect. 1, the ’’background” amplitudes and behave 

like l/w' in the low energy limit, while is proportional to to'. Therefore, 

the expansion of the spin-averaged matrix element reads 

-| ^BH~\-Born ^BH~\-Born 

- = + + + + (14) 

spins 


where the coefficients C^H+Bom^ fjBH+Bom (jBH+Bom entirely de¬ 
termined by the background. The coefficient Cq°^ results from the interference 
of the singular part of the background term with the leading term of 
and contains the information on the GPs. It can be expressed in terms of four 
structure functions lil. 


Cq°’' = K2{e,q)^vi{ePLLiq) - PTT{q)) -f U2 V2e(l-f e)PLT(g) 
+'i'3\/2e(l -I- ^jP^riq) I 


(15) 


with kinematical factors Vi depending on g, 9, (j) and e. The structure functions 
contain the GPs in comhmation with the elastic form factors of the nucleon. 
Using the results of Ref. Ell, the structure functions Plt and are mutually 
dependent via the relation 


Phriq) 


—p^PLTiq) ■ 

Wo 


(16) 
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As a consequence, the coefficient C(f°* contains only three independent struc¬ 
ture functions, thus leading to information on three of the six independent 
GPs. 

In order to disentangle the remaining polarizabilities one has to resortAp 
double polarization observables. As has been suggested by VanderhaeghenEj, 
the reaction p(e, e’p')^ can be used to extract the remaining three polarizabil¬ 
ities. Instead of measuring recoil polarization, the same information can be 
gained by the reaction p{e,e'p' )j. We focus the attention on the observables 
for beam-target polarization in order to give some guidance for a possible ex¬ 
periment with the BLAST facility at MIT-Bates. To this end we consider the 
target (double) asymmetries 


= 7 -- ^7 -^^ , (17) 

where h is the electron helicity and rrii (i = x, y, z) is the spin-projection of 
the nucleon target with respect to the coordinate system 


q' — cos dq ^ q X q' 
sin ^ sin d 


(18) 


As in the unpolarized case, one can perform a low energy expansion for the 
numerator of the asymmetries leading to expansion coefficients In 

terms of structure functions these coefficients read 


^o{to) =4A:2(e,g) ?;i V^W^7TT°^(g) -f 


+V 3 \/2e(l-e)P^^y’°^ (q) 


LT 

= 4if2(e, q) [< v'2e(l - (g) + (?) 

. 0 ) 


+V 3 \/l — Pr^rp ’ \q) + \/2e(l — e)Ppp 


(19) 


( 20 ) 


A similar decomposition of contains the same structure functions as 

^o{xO) ^9- (10), however with new kinematical factors wf. Note that the 


kinematical factors for target polarization in z-direction are the same as in 
the unpolarized case. Altogether one obtains three new independent structure 
functions, extraction of requires an 

out-of-plane measurement since uf —> 0 for (j) 0°. This situation is quite 

analogous to the case of double polarization asymmetries with recoil polariza- 
tionllo, whose expansion coefficients contain the same kinematical factors as 
for target polarization, while the corresponding structure functions may have 
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Figure 3: Background cross section for the Bates kinematics and the MAMI kinematics as 
function of the scattering angle 9. Dashed line: VCS Born contribution, dotted line: Bethe- 
Heitler contribution, full line: complete background cross section. 


a somewhat different appearance. The explicit expr^ions for the coefficients 
^o{fo) presented in an upcoming publication CJ. 

3.2 Results and Discussion 

In the following numerical studies we concentrate on two different kinematics: 
(I) the Bates kinematics with q = 240 MeV, e = 0.9 and lo' = 100 MeV and (II) 
the MAMI kinematics with q = 600 MeV, e = 0.62 and lo' = 111 MeV. Bot i 
sets of vari^les have been explored in the unpolarized experiments at Bates □ 
and MAMIB. The main difference beteween the two kinematics is in the values 
of q and e, while the difference in to' will be neglected in the discussion below. 
In our analysis, the observables are always calculated as function of the scatter¬ 
ing angle 6. In addition, we distinguish between in-plane {(p = 0) and the 90° 
degree out-of-plane kinematics for the Bates kinematics, while in the case of 
the MAMI kinematics the in-plane and 60° degree out-of-plane cases have been 
investigated. Results for the background cross section are shown in Fig. H, with 
the separate contributions from the BH and the VCS Born amplitude. Bor the 
elastic form factors of the proton we used the parametrization of Ref. E3. In 
the case of in-plane kinematics, the angular dependence of the BH cross sec- 
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Figure 4: Influence of the GPs on the unpolarized cross section of the reaction p{e, e'p' )-f 
for Bates and MAMI kinematics. Full line: LSM result, dash-dotted line: HBChPT result. 
The quantity a — crsH+Bom represents the contribution of the coefficient Cq°^ of Eq. ( p^ 
to the cross section. 


tion is characterized by two peaks occurring for the real photons emitted along 
the direction of the initial or final electrons. In these regions, the BH term 
gives the overwhelming contribution to the background cross section, while it 
becomes much smaller at negative values of 9. For out-of-plane configurations 
the two peaks disappear, and in the case of the Bates kinematics the angular 
distribution of the background term is almost symmetrical around 0 = 0. In 
order to increase the sensitivity to the GPs, the region of negative scattering 
angles corresponding to photon emission in the half-plane opposite to the out¬ 
going electron has to be explored. In comparison to the MAMI kinematics, the 
background cross section is one order of magnitude higher in the Bates case, 
whereas the 9 dependence of the background is similar for both kinematics. 

For the LSM and HBChPT the relative influence of the GPs on the unpo¬ 
larized cross section has been plotted in Fig. As is to be expected, HBChPT 
generally predicts larger effects than the LSM. In both calculations the influ¬ 
ence of the GPs is negative for a large region of 9, leading to a reduction of the 
total cross section. For the unpolarized case the anomaly gives no contribution 
as long as only the leading order in lo' is taken into account. 

For the Bates kinematics the GPs change the cross section by about 20%, the 
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effect being almost completely given by a{q). The influence of the GPs for the 
MAMI kinematics is of the order of 15%. In this case the LSM predicts that 
only 60% of the signal are due to the electric polarizability. Since one aims at 
cross section measurements with a relative error of about 2%, the calculated 
effects should be detected. However, there are two problems which render 
the measurement of the GPs more difficult: (I) the effect of the GPs strongly 
increases with uj'. At cu' = 70MeV, e.g., the effects shown in Fig. ^ reduce 
by roughly 50%. For higher values of uj', where the signal is larger, higher 
order terms of the residual amplitude can no longer be neglected (see also the 
corresponding discussion in sect. 2.1). In this case the challenging task is the 
reliable extraction of the leading order contribution in which one is interested. 
(II) Independent of the value of u' there are large radiative corrections c3. 
These corrections result in a 20% contribution to the cross section which is of 
the same order of magnitude as the influence of the GPs, thus complicating 
the separation of background and GP effects. 

In Fig. H we show the GP contribution to the target asymmetries and 
for both LSM and HBChPT, without taking account of the anomaly con¬ 
tribution. Obviously, larger effects can be expected in the case of . While 
in the unpolarized cross section the GP effect was larger for Bates kinematics, 
we now observe the opposite situation for the asymmetries. Here the maxi¬ 
mum signature is about 10% for the MAMI kinematics but only about 3% for 
the Bates kinematics. The different behaviour between the unpolarized and 
the polarized case has the following explanation: In Fig. ^ we have seen that 
the background, which determines the denominator of the asymmetries, is one 
order of magnitude larger for the Bates kinematics. On the other side, the 
influence of the GPs in the numerator of the asymmetries increases much less 
when going from Bates to MAMI kinematics, while the influence of the GPs on 
the unpolarized cross section also strongly increases by changing the kinemat¬ 
ics, in particular due to the strong q dependence of thejSlectric polarizability. 
More details of this discussion will be presented in Ref. E£l. 

In Fig. ^ we display the anomaly contribution to the asymmetries as predicted 
by the LSM. Since the influence of the anomaly is strongly decreasing with in¬ 
creasing q, the largest effect shows up for the Bates kinematics. In this case the 
shape of the asymmetry is completely changed by the anomaly. Though 
for the MAMI kinematics the anomaly is less dominant, its influence is still 
quite sizeable. 

We conclude that a full determination of the GPs will require double po¬ 
larization experiments, polarized electrons and target or recoil polarization, 
in addition to the unpolarized experiments. Given the projected accuracy of 
the experiments, one should be able to measure the GPs in spite of a large 
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Figure 5: Influence of the GPs on the target asymmetries for Bates and MAMI kinematics. 
Full line: LSM result, dash-dotted line: HBChPT result. The quantities 

and , d represent the contribution of the coefficients and in 

_ _I BiUsBorn ^ 0(z,0) 0(a:,0) 

Eqs. ( |l£| ) and (|20[), respectively. 
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Figure 6: Target asymmetries with and without the anomaly for Bates and MAMI kinematics 
as obtained in the LSM. Full line: anomaly included, dashed line: without anomaly. 











background. Since the momentum dependence of the GPs is predicted quite 
differently by various models, such precision experiments will be invaluable to 
restrict the model parameters, and, quite generally, as further benchmarks of 
the nucleon’s structure. 
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